Superconducting fluctuations in the Luther-Emery liquid 



Edmond Orignac 

Laboratoire de Physique Theorique de I'Ecole Normale Superieure CNRS-UMR8549, F-75231 Paris Cedex 05, France 

Didicr Poilblanc 

Laboratoire de Physique Theorique CNRS-FRE2603, 
Universite Paul Sabatier, F-31062 Toulouse, France 
(Dated: February 2, 2008) 

The single-particle superconducting Green's functions of a Luther-Emery liquid is computed by 
bosonization techniques. Using a formulation introduced by Poilblanc and Scalapino [Phys. Rev. B 
66, 052513 (2002)], an asymptotic expression of the superconducting gap is deduced in the long 
wavelength and small frequency limit. Due to superconducting phase fluctuations, the gap exhibits 
as a function of size L a (l/L) 1//2Ap power-law decay as well as an interesting singularity at the 
spectral gap energy. Similarities and differences with the 2-leg t-J ladder are outlined. 
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Superconductivity in cuprate materials has sparked an 
important theoretical effort to investigate superconduct- 
ing order in strongly correlated materials. A long-ranged 
superconducting ground state (GS) is characterized by 
a non-zero superconducting Gorkov's off-diagonal one- 
electron Matsubara Green functional defined as 



F(q,r) = -<T C _ q ,_ CT (r/2)c q , CT (-T/2)) 



(1) 



As proposed recently Q, basing on the Nambu- 
Eliashberg theory 0, the Matsubara frequency- 
dependent gap function Agrj(q, ui n ) can be directly 
computed from the diagonal and off-diagonal Green 
functions G(q, u n ) and F(q,u) n ) as 



Asc(q,w„ 



2iw„F(q, Lu r , 



G(q,w n ) - G(q.,-u> n ) 
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Such an expression has been used in numerical computa- 
tions for the two-dimensional t-J model0] and for two-leg 
t-J ladders 0. A one dimensional system such as the two- 
leg t-J ladder is actually far remote from a normal Fermi 
liquid or a conventional superconductor 0,0, 0. The t- 
J ladder belongs to the the class of Luther-Emery (LE) 
liquids 7], which are characterized by spin/total charge 
separation, with a spectral gap in spin excitation and 
gapless charge excitations 0. The gapless charge excita- 
tions lead to a perfect metallic conductivity [9| without 
any long range superconducting order 0. The spectral 
gap in the spin excitationsgenerates a gap in the one- 
electron spectral functions|llill2|. In an infinite system, 
the absence of long range superconducting order implies 
that the correlation function has to vanish. However, 
in a finite system, the function Q does not vanish and 
the quantity defined in may provide useful informa- 
tion on the pairing processes 3]. As a complementary 
study to numerical studies, it is interesting to investigate 
analytically the behavior of A$c (q, w ) in a finite size LE 
liquid. A complication in the case of the two-leg ladder 



is that interband (but not total) charge excitations con- 
tribute to the formation of the spin-gap Q. A simpler 
example of a LE liquid is afforded by the U < Hub- 
bard chain IT3. IT^ or the related spin-gap phase of the 
t-J chain 15] in which the gap is developed purely in 
the spin mode. Moreover, this model of a LE liquid is 
integrable both on the lattice [l4| and in the continuum 
limit ^(|. Progress in the theory of integrable systems 
has made it possible to calculate correlation functions 
using the form factor expansion^]- In the case of gap- 
ful systems, it has been shown that the first few terms of 
the expansion usually lead to good approximations of the 
physical quantities[18|. In the present paper, we will cal- 
culate Asc in the case of an attractive one-dimensional 
Hubbard chain. The Hamiltonian reads: 
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H. c.) + U 'Y^ n p,t n p,l- ( 3 ) 



Although the Hubbard model in one-dimension is in- 
tegrable, its form factors have proved very difficult to 
obtain[lJ|- Since we are mostly interested in the low- 
energy and long-wavelength, it is more convenient to 
take the continuum limit. The continuum Hamiltonian 
H = H„ + H a reads: 
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H a = 



dx 
o 2^ 
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where </> p and 



are the bosonic fields describing respec- 
tively the spin and charge excitations, n p and Tl a their 
conjugate variables, a is a short distance cut-off (~ lat- 
tice spacing) and L ~ Na is the system size. Note that 
we assume here a unique velocity v for spin and charge 
excitations. The spin excitations are described by the 
sine- Gordon model. The term cos \f&4><y 1S marginally 



relevant and opens a spin gap M. . The excitations above 
the ground state are known |l6| to consist only of mas- 
sive solitons (spinons) carrying the spin ±1/2 and charge 
zero, in agreement with the exact solution of the lat- 
tice Hubbard model 0]. ^ is important to note that 
with the usual definition of the spin gap A a as the gap 
between the ground state and the lowest triplet excited 
state, A ff = 2M.. The fermion operators can be decom- 
posed as c P: a — \faY^ r e%kFrx 4 } r,a{x) (x — pa, r = ± or 
±1) and the fermionic fields ip r>a can be expressed as 0: 



ip r ,A x ) 
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where 9 P and Q a are the dual fields. Fermion propaga- 
tors thus factorize into the product of a holon propagator 
and a spinon propagator. Since the charge Hamiltonian 
is quadratic, the calculation if the holon propagator is 
straightforward. For the spinon propagator, one needs 
to know the expression of the form factors of soliton- 
creating operators e ± l{ fi°- r <t>v) m the sine-Gordon model. 
These have been obtained recently in infinite volume [20j 
and these results have been applied to the calculation of 
the diagonal Green's functions |l2| confirming the ansatz 



of Ref. 11]. In our case, we need form factors in a system 
of finite size. However, if the system size L is much larger 
than the correlation length £ = v/M, since the dominant 
contribution to Ascil-, comes from integration on the 
region x, vr <C £, it is a good approximation to replace 
finite volume form factors by infinite volume ones. Keep- 
ing this in mind, we define the diagonal and off-diagonal 
Green's functions as: 

G(x,r)=- (T T W,. (x, t)^ (0, 0)) , (6) 

F(x,T) = -<T T Vr,a(iB,T)^_ P ,_ ff (0 ) 0)>. (7) 

Following ^3 , the asymptotic (long-distance) behaviors 
of the diagonal and off-diagonal Green's functions is then 
obtained as, 
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where p = ^x 2 + {vt) 2 , V = \{K p + l/K p - 2), 77' = 
\(K n — l/Kp), Z\ is a dimensionless constant calculated 



m l^llj and M the single particle gap. These expressions 
show that G and F have a similar exponential decay 
at large distances above the characteristic length-scale 
£ = v/M.. More rapidly decaying terms like e~ 3Mp / v 
and higher that involve more than one spinon in the in- 
termediate state have been neglected. Note that the SC 

1 

Green's function decays with system size like as (l/£) 2Kp 
due to SC phase fluctuations. However, apart from this 
overall power-law scaling factor, we expect the SC gap 



FIG. 1: Superconducting gap (normalized to its Q = value) 
vs Q£. 



Age to remain finite and bear interesting q and u) de- 
pendence. Using the Fourier transforms of the above ex- 
pressions (0, the momentum and Matsubara- frequency 
dependence of the gap is given by 



A sc ( q ,u; n )=CM(^y 



2 F x (f -r/,f-7/;l; y) 



2j Fi(1 - 77 
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(9) 

where y= -(Q£) 2 , Q 2 = (q - k F ) 2 + (uj n ) 2 /v 2 , 2 F 1 is the 
hypergeometric function |2l| and C = T(| — 2?/)/r(2 — 
2rj). The approximate expression 10 although a long- 
wavelength limit still applies up to momentum Q ~ 27r/£. 
Note that (Qa) 2 terms are dropped since £ ^> a. Asc 
is maximum (minimum) at Q — (i.e. at the Fermi mo- 
mentum) for K p < 1.57 (K p > 1.57) and shows only a 
very moderate dependence in Q as shown in Fig. ^ F° r 
Q -> 0, A sc (<Z,Wn) =s A SC (0,0)(1-C'(Q£) 2 ). The con- 
stant C = h \\K 2 - AK a - 13 



lK - 



1 



decreases 



for increasing SC correlations i.e. for increasing K p , e.g. 
C ~ 0.674, 0.188 and 0.016 for K p = 0.5 (charge den- 
sity wave regime), 1 and 1.5 (superconducting regime), 
respectively. 

The real-frequency expression of the gap is obtained 
by the analytic continuation (w„) 2 — > — (ui + it) 2 in 
Eq. ©. For |cj| > M, the hypergeometric functions 
have branch cuts, leading to a nonzero imaginary part 
in Asc(0,cj) as can be seen on Fig. Hfb). Above the 
threshold, the gap function presents a singular behavior: 

A(0,w) ~ A(0,0)(1 - {ul ^ )2 y^ ■ The divergence at 
u> = M is in fact cut-off once \ui — M\ < 1/L since we 
are really dealing with a system of finite size, so that 
Ascfew — > M) = M. The full w-dependence of Asc 
plotted in Fig. Ola) indeed reveals a strong singularity 
at lj = M increasingly pronounced as system size is in- 
creased. 

In the absence of data for a single attractive Hub- 
bard chain or a single t — J chain, we compare our re- 
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This implies a second threshold in A(q, ui n ) at the en- 
ergy oj = M(l + V2) besides the threshold at energy 
a; = M.. Whether this could be related to some higher 
energy features seen in numerics Q needs further clar- 
ifications. This will be discussed in more details in a 
separate publication. 

In conclusion, we have computed the fluctuating SC 
gap of the LE chain. A simple form is obtained with a fac- 
torization into a power-law factor accounting for SC sup- 
pression due to quantum phase fluctuations multiplied 
by a function containing the dynamics of the pairing in- 
teraction. We point out some differences with the case of 
the 2-leg t-J ladder and suggest that the gapped sectors 
of the latter could be better described by the SO (6) GN 
model. 



co/M 



FIG. 2: Real (a)) and imaginary (b) parts of the supercon- 
ducting gap vs oj obtained by bosonization results for a ID 
LE liquid (energies are in units of the spectral (solitonic) gap 
A4) for K p = 1.3 and q = Icf. Data for two lengths are shown. 



suits to numerical calculations on the t-J ladder model 
at 1/8-doping and J = 0.4 3j. The numerical calcula- 
tions on 2 x 12 ladder show that there are two spinon 
gaps, M(q y = 0) = 0.08i and M(q y = %) = 0.12f such 
that F(q y = 0, tt), G(q y = 0, tt) and A(q y = 0,ir) de- 
velop an imaginary part for oj > A(q y ). The compari- 
son of these s pin on gaps with the actual spin gap of a 
2 x 24 ladder |24| (defined as the gap between the low- 
est triplet excitation and the singlet ground state) gives 
A a {q y = 0) = O.lli ~ 2M(q y = 0) and A a (q y = tt) = 
0.17* ~ M(q y = 7r) + M(q y = 0). The discrepancies 
could result from a spinon-spinon attraction or from the 
overestimation of the spinon gaps in the 2 x 12 ladder. 
It is also important to note that in the correlation 
length is of the order of magnitude of the system size. 
We note that no sharp peak is present in the imaginary 
part at the threshold, in contrast to the prediction of 
0, but the prediction of a rather constant behavior of 
A below the threshold is in agreement with 10. 

A more detailed comparison between analytic and 
numerical result is possible. In the case of the lad- 
der system, away from half-filling, the gapped modes 
is expected to present an approximate 50(6) ~ St/ (4) 
symmetry |22l l2fi| . This allows a description of the 
gapped modes by the SO(6) Gross-Neveu (GN) model^ 
and a form factor calculation of the superconducting gap 
A(q, uj n ) along the lines of the present paper [2(J. The 
novelty in the case of the SO(6) GN model is that on top 
of the spinon excitations of mass Ai (known as kinks in 
the literature on the GN model), there are also massive 
fermion excitations (bound states) with a mass A4\^2. 



[» 

[9 
[10 

[11 
[12 

[13; 

[14 
[IS 

[ie; 

[17 

[is; 



J. R. Schrieffer, in Theory of Superconductivity (Addison- 
Wesley, Reading, MA, 1964), Chap. 7. 
D. Poilblanc and D. J. Scalapino, Phys. Rev. B 66, 
052513 (2002). 

D. Poilblanc, D. J. Scalapino and S. Capponi (unpub- 
lished) . 

E. Dagotto and T. M. Rice, Science 271, 618 (1996), and 
references therein. 

H. J. Schulz, in Correlated fermions and transport in 

mesoscopic systems, edited by T. Martin, G. Montam- 

baux, and J. Tran Thanh Van (Editions frontieres, Gif 

sur Yvette, France, 1996), p. 81. 

E. Dagotto, Rep. Prog. Phys. 62, 1525 (1999). 

H. J. Schulz, in Mesoscopic quantum physics, Les 

Houches LXI, edited by E. Akkermans, G. Montambaux, 

J. L. Pichard, and J. Zinn- Justin (Elsevier, Amsterdam, 

1995), p. 533. 

For numerical studies see e.g. C. A. Hayward et 
al, Phys. Rev. Lett. 75, 926 (1995); D. Poilblanc, 
D.J. Scalapino, and W. Hanke, Phys. Rev. B 52, 6796 
(1995). 

H. J. Mikeska and H. Schmidt, J. Low Temp. Phys 2, 371 
(1970). 

N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 
(1967). 

J. Voit, Eur. Phys. J. B 5, 505 (1998) 

A. M. Tsvelik and F. H. L. Essler, |cond-m at/0205294 
(unpublished). 

A. Luther and V. J. Emery, Phys. Rev. Lett. 33, 589 
(1974). 

N. Andrei, in Low- Dimensional Quantum Field Theories 
For Condensed Matter Physicists, edited by S. Lundqvist, 
G. Morandi, and L. Yu (World Scientific, Singapore, 
1993), and references therein. 

M. Nakamura, K. Nomura, and A. Kitazawa, Phys. Rev. 
Lett. 79, 3214 (1997) , Eond-mat/9708204] 
H. Lc - 



1698 (1979). 

M. Karowski and P. Wiesz, Nucl. Phys. B 139, 455 
(1978). 

D. Controzzi, F. H. L. Essler, and A. M. Tsvelik, in 
New Theoretical approaches to strongly correlated sys- 



4 



terns, Vol. 23 of NATO Science Series II. Mathematics, 
Physics and Chemistry, edited by A. M. Tsvelik (Kluwer 

Academic Publishers, Dordrecht, 2001), p. 25. 

[19] F. H. Essler and V. E. Korepin, cond-mat/9808018 (un- 
published). 

[20] S. Lukyanov and A. B. Zamolodchikov, Nucl. Phys. B 

607, 437 (2001). 
[21] M. Abramowitz and I. Stegun, Handbook of mathematical 

functions (Dover, New York, 1972). 
[22] H. Schulz, cond-mat/9808167 (unpublished). 
[23] D. J. Gross and A. Neveu, Phys. Rev. D 10, 3235 (1974). 



[24] Data obtained using an approximate Contractor Renor- 
malisation (CORE) method; see S. Capponi and D. Poil- 
blanc, Phys. Rev. B 66, 180503 (2002) and references 
therein. 

[25] Numerical investigations of boundstates are given in 

D. Poilblanc, O. Chiappa, S. White and D. J. Scalapino, 
Phys. Rev. B 62, R14633 (2000). 

[26] A calculation of the spectral functions in the ladder using 
the SO(6) GN model has been reported in C. Ahn and 

E. Paik, J. Korean Phys. Soc. 39, 965 (2001). 



